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INTRODUCTION: 


Many  workers  have  applied  column  chromatography  to  groundwater 
contaminant  transport  problems  and  several  have  missed  the  conclusion  of 
Bear  (1979),  who  states  that,  "because  of  both  technical  and  scaling 
difficulties,  the  sand  box  model  cannot  be  used  as  a  tool  for  solving 
boundary  value  problems  of  practical  interest.  The  column  filled  with  a 
porous  medium,  however,  is  a  very  important  laboratory  research  tool  suited 
for  observing  phenomena  (also  of  solute  transport  with  or  without  solid- 
liquid  interaction)  and  verifying  approximate  solutions".  Investigations 
reported  here  on  sorption  coefficient  measurement  of  single  sorbates  and 
sorbate  mixtures  on  aquifer  material  columns  were  done  in  accordance  with 
this  caveat. 

In  this  work,  chromatography  columns  were  packed  with  aquifer  material 
and  sorption  coefficients  were  determined  by  frontal  analysis.  This 
Chromatographic  measurement  of  equilibrium  sorption  was  developed  by  James 
and  Phillips  (1954)  among  others,  elaborated  by  Jacobson  et  al  (1984),  and 
applied  to  soils  and  aquifer  materials  by  eg.  Bilkert  and  Rao  (1985), 
Schwarzenbach  and  Westall  (1981),  Goerlltz  (1984),  Kay  and  Elrick  (1967), 
Russell  and  McDuffie  (1986).  It  has  also  been  used  on  synthetic  sorbates 
(e.g.,  Sharma  and  Fort  1973).  Chromatographic  measurements  of  equilibrium 
sorption  are  generally  referred  to  as  dynamic,  in  contrast  to  conventional 
batch  (static)  equilibrium  sorption  measurements.  Under  ideal  conditions, 
equilibrium  sorption  isotherms  for  a  given  solvent-sorbate-sorbent  system 
should  be  independent  of  whether  the  isotherm  is  determined  by  a  dynamic  or 
static  method.  The  thermodynamics  of  sorption  from  simple  solution  has  been 
presented  by  Everett  and  Podoll  (1979),  who  also  consider  isotherm 
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measurements  by  dynamic  and  static  methods.  In  the  present  work, 
equilibrium  sorption  data  on  water-nonionic  organic  compound-aquifer 
material  systems  are  determined  by  both  dynamic  and  static  methods.  Results 
are  compared  and  reasons  for  differences  between  the  methods  in  this 
application  to  aquifer  materials  are  discussed.  Conditions  for  ideal  liquid 
column  chromatography  with  frontal  analysis  are  given  by  Snyder  and  Kirkland 
(1979).  Among  these  are:  constant  flow  velocity,  uniformity  of  packing 
particle  size,  homogeneous  packing,  absence  of  column  wall  and  end  effects, 
complete  retention  of  packing  by  a  packing  retainer,  flow  velocity  low 
enough  to  attain  local  equilibrium  of  solute  between  solvent  and  packing  but 
high  enough  so  dispersion  does  not  override  sorption  effects.  Some  of  these 
conditions  cannot  be  met  with  aquifer  material  packing,  so  they  must  be 
approached  as  closely  as  possible. 

Transport  of  solute  through  a  column  is  well  known  for  the  case  of 
solute  sorption  with  a  linear  isotherm.  The  sorbates  used  here,  1- 
methyl naphthalene  and  1  ,2-dichlorobenzene  have  linear  isotherms  on  the 
aquife”  materials  selected.  The  operation  of  a  finite  length  (L)  ideal 
column  Is  described  following  Bear  (1979)  for  the  case  of  linear  sorption  in 
terms  of  the  one-dimensional  mass  balance  equation. 
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where  C  is  the  solute  concentration  in  the  liquid  phase,  D  is  the  axial 


dispersion  coefficient,  Uo  is  the  liquid  flow  velocity,  k  is  the  sorption 
coefficient,  Z  is  distance  (axial)  and  t  is  time,  is  packing  bulk  denisty 
and  e  is  porosity. 

For  frontal  chromatography  the  initial  and  boundary  conditions  are: 

at  t  <  0  for  z  >  0,  C  =  0,  initial  condition; 

9c 

for  t  >  0  for  z  =  0,  (C^-C)  =  -D  ,  inlet  boundary  condition; 

3C 

z  =  L,  _  =  0,  outlet  boundary  condition. 

3z 

The  inlet  boundary  condition  corrects  the  concentration  gradient  for 
dispel  Sion,  while  the  outlet  boundary  condition  implies  flow  into  a  gas 
Continuum.  Bastian  and  Lapidus  (1956)  have  provided  an  analytic  solution 
for  this  situation  (i.e.,  a  functional  relation  for  C(Z,t)/C^,  where  is 
the  column  feed  concentration. 

The  column  output  response  is  affected  by  dispersion  of  the  solute 

D 

indicated  bv  the  term  _  _  and  by  advection  of  the  solute 
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indicated  by  the  term  ® 
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Both  terms  contain  the  equilibrium 
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sorption  coefficient.  For  constant  C,  Uq,  D,  and  k,  the  breakthrough  curve 
for  C/C^  versus  time  at  the  outlet  is  sigmoid  as  shown  in  Fig.  1.  Frontal 
elution  occurs  at  time  t^  for  an  unretained  solute  and  t^  for  a  retained 
solute . 

The  difference  in  area  between  the  two  curves  in  Fig.  1  is  a  measure  of 
the  retained  solute  sorption  coefficient  k.  Letting  be  the  concentration 
of  unretained  solute  and  be  the  concentration  of  the  retained  solute, 
with  the  subscript  (0)  indicating  input  concentration,  the  elution  time  t^ 

1  r® 

equals  _  Jo  C^dt,  and  the  retention  volume  equals  Uotj. 

(0) 

A  similar  expression  determines  V^,  the  retention  volume  of  the  sorbed 
solute.  A  point  on  the  isotherm  corresponding  to  sorbate  equilibrium 
concentration  is  calculated  as  given  by  Jacobson  et  al  (1984),  namely, 


q,  _  V,  -V 


Where  q^  =  equilibrium  concentration  of  sorbate  on  sorbent  at  and  Vsp 

is  the  volume  of  sorbent  in  the  column.  Vsp  is  usually  determined 
grav imetr i cal ly .  Vj  is  often  called  the  column  dead  volume.  Units  used 
here  are  in  mmol  and  q^  in  mmol/ml  of  sorbent. 

There  are  several  practical  difficulties  that  must  be  confronted.  A 
complete  isotherm  by  frontal  analysis  requires  repeat  column  runs  at  several 
^2(n)  establish  isotherm  points.  This  is  a  slow  process  because 

u^  must  be  low  enough  to  achieve  local  equilibrium.  No  retainer  could  hold 
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An  example  of  normalized  breakthrough  curves  for  retained  and  unretained 
solutes . 


up  the  aquifer  fines  without  acting  as  a  sorbent  and  producing  excessive 
backpressure  at  practical  flow  rates.  There  is  no  way  to  homogeneously  pack 
aquifer  materials  in  a  column  because  of  non-uniform  particle  size. 
Channeling  may  occur  and  not  be  detected.  Plug  flow  may  not  occur  due  to 
Column  wall  and  end  effects.  Flow  may  be  too  fast  for  local  equilibrium 
attainment.  There  is  no  completely  unretained  solute  to  measure  ,  so  the 
least  unretained  solute  must  be  used.  Sorption  by  the  chromatographic 
apparatus  may  cause  error  in  the  calculated  sorption  coefficient. 

Previous  researchers  using  soils  and  aquifer  materials  did  not  state 
column  packing  retainer  characteristics.  Concern  developed  here  when 
aquifer  material  fines  were  noted  to  pass  through  the  2  micron  stainless 
steel  frits  used  in  trial  column  loading  and  flushing  tests  (and  throughout 
the  later  experimentation).  Goerlitz  (198i<),  Wilson  et  al  (1981),  Russell 
and  McDuffie  (1986)  and  Nkedi-Kizza  et  al  (1983)  reported  experiments  on 
columns  packed  with  soil  or  aquifer  material,  but  did  not  adequately  specify 
the  retainer  material  and  pore  size.  There  is  an  apparent  lack  of  concern 
with  particle  retention,  but  the  effects  of  particle  losses  on  sorption 
coefficients  cannot  now  be  estimated.  The  aquifer  materials  used  here  were 
selected  from  a  number  of  geological  settings  and  all  contained  considerable 
fines.  It  is  assumed  that  the  loss  of  fines  will  be  a  problem  in  column 
work  with  most  aquifer  materials,  and  this  is  now  being  investigated. 

Materials 

CaCl22H20  (Baker  Reagent  Grade)  was  used  to  prepare  CaClj  solutions. 

Reagent  grade  1 -methylnaphthalene  (ICH3N)  and  1,2  dichlorobenzene 
(ODCB)  of  99  percent  purity  were  used  as  obtained  from  the  suppliers. 
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Aldrich  Chemical  Co.,  Inc.,  Milwaukee,  WI ,  and  Fluka  Chemical  Co., 

Ronkonkoma,  NY,  respectively.  Water  was  MiliQ  (Milipore  Corp.)  treated 
after  distillation  in  glass.  Solvents  used  for  cleaning  apparatus  and 
extraction  were  Burdick  and  Jackson  distilled-in-glass  CHjCl,  acetone  and 
hexane  from  American  Scientific  Products  Co.  JP-4  jet  fuel  was  supplied  by 
the  Defense  Logistics  Agency  Fuel  Depot  at  Lynn  Haven,  FL . 

Aquifer  materials  were  obtained  from  several  sources.  Lula  material 
was  collected  at  a  depth  of  5m  from  the  Johnson  Ranch  near  Lula,  OK,  and 
provided  by  Robert  S.  Kerr  Environmental  Research  Laboratory  (EPA).  Ada,  OK. 
The  USAF  provided  aquifer  materials  from  its  bases  (which  are  used  as  the 

sample  names).  Tinker  material  was  collected  near  Oklahoma  City,  OK,  from  a 

depth  of  4.9  -  8.7m.  Barksdale  aquifer  material  was  collected  near 
Shreveport,  LA,  from  a  depth  of  4.3  -  13.9m.  Blytheville  aquifer  material 

Was  collected  near  Blytheville,  AR,  from  a  depth  of  8.8  -  10.0m.  Carswell 

material  was  collected  near  Fort  Worth,  TX,  from  a  depth  of  1 .3  -  6.4m. 
Canadian  Borden  site  material  was  a  subsample  of  the  bulk  sample  provided  by 
Mackay  et  al  (1986).  Characteristic  properties  and  mineralogi cal 
description  of  these  samples  are  given  in  detail  by  Stauffer  (1987)  and 
selected  values  are  presented  in  Table  1. 

Table  1.  Aquifer  Material  Properties. 


Aquif  er 
Material 

TOC  (%) 

CEC 

(mcq/1 OOg) 

Surf  ace 

Area 

(m,/g) 

Sand 

i%) 

Clay 

i%) 

Silt 

{%) 

Iron 

{%) 

T inker 

0.010 

20.9 

9.2 

40.8 

4.7 

53.1 

1.4 

Carswell 

0.027 

8.0 

9.5 

74.8 

NR 

18.2 

0.9 

Barksdale 

0.105 

32.4 

7.5 

52.3 

NR 

41  .5 

1  .0 

Blytheville 

0.156 

18.1 

8.0 

67.8 

3.5 

27.3 

0.7 

8 


Border. 

0.015 

NA 

0.3 

96.0 

2.0 

2.0 

1  .5 

Lula 

0.020 

6.3 

11.8 

91  .0 

s.ii 

5.6 

2.9 

NA  =  not  analyzed. 

nr  =  not  measurable  in  sieve. 

Experimental  Methods. 

The  ohromatographic  apparatus  was  constructed  of  stainless  steel  and 
glass  to  minimize  sorption  problems,  and  is  shown  in  Fig.  2.  All  columns 
used  were  25  cm  long  stainless  steel  liquid  chromatography  blanks  with 
stainless  steel  packing  retainers  supplied  by  Supelco,  Bellfonte,  PA. 

Column  retainer  frits  were  2  micron  pore  size.  The  switching  valve  was  a 
low  dead  volume  manual,  6  port,  Nitronic-60  alloy  Valeo  unit  with  1/16"  zero 
volume  fittings  supplied  by  Supelco.  All  tubing  and  connections  were 
stainless  steel.  Tubing  ID  was  0.010"  ID  used  in  minimum  lengths. 

Detectors  were  an  ISCO  variable  wavelength  uv  model  and  a  Perkin-Elmer 
r’efractive  index  model  LC-25.  A  Gilson  (Middleton,  WI)  model  FC-80N 
fraction  collector  was  used  to  collect  portions  of  column  effluent  for 
analysis  of  Cl  by  titration  with  silver  nitrate. 

Chromatographic  determination  of  a  sorption  isotherm  point  required 
about  100  hours.  Detector  signals  were  recorded  on  a  Hewlett-Packard 
3357  laboratory  data  system.  Files  of  5*10  min  duration  were  taken, 
transferred  to  the  Hewlett-Packard  RPN  program  and  merged  to  make  a  complete 
file  for  a  chromatography  run  Frontal  analysis  integrations  were  done  on 
the  merged  files  by  commands  in  the  RPf^  system. 

Frontal  analysis  was  done  after  preparation  of  a  saturated  solution  of 
the  organic  sorbate  in  water.  Saturation  vessels  and  techniques  used  in 
sorbate  solution  preparation  are  described  in  Burris  and  MacIntyre  (1985). 
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Solutions  were  prepared  in  0.01 5M  CaCl^-  Experiments  were  conducted  at  a 
laboratory  temperature  of  22°Cj^2°C.  The  saturated  solution  was  then 
transferred  to  a  ^lass  mixing  chamber  by  positive  pressure  as  shown  in  Fig. 
3.  Dilution  water  0.015M  in  GaCl^  was  added  to  make  up  the  desired  solution 
concentration.  The  contents  of  this  mixing  chamber  were  transferred 
immediately  to  the  HPLC  pump  2  chamber.  Occasional  checks  of  the  sorbate 
concentration  were  made  on  20  mL  portions  taken  from  the  pump  delivery  as 
described  by  Burris  and  MacIntyre  (1986).  Concentrations  of  sorbates  are 
reported  as  fractions  of  their  saturation  concentration  value,  and  values 
above  0.75  times  saturation  concentration  were  avoided  due  to  potential 
problems  with  sorption  on  the  apparatus  or  nonlinearity  in  sorbate 
isotherms.  Pump  1  was  filled  with  0.01M  CaCl^.  Both  pumps  were  set  to  the 
same  flow  rate  and  the  column  was  flushed  till  the  solute  was  free  of 
visible  particles  and  the  uv  detector  output  was  steady.  The  valve  was  then 
rotated,  producing  a  step  increase  from  0.01M  to  0.015M  in  CaC^ 
concentration  and  from  zero  to  the  specified  organic  sorbate  concentration. 
This  generated  fronts  for  analysis.  The  CaCl2  breakthrough  was  detected  by 
refractive  index,  and  later  by  fraction  collection  and  AgNO^,  titration  of 
each  fraction.  Breakthrough  of  organics  was  detected  by  uv  absorbance  at 
28l  nm  for  1-methylnaphthalene  and  268  nm  for  1,2-dichlorobenzene. 
Confirmation  of  organic  breakthrough  was  determined  by  manual  collection  of 
samples  at  regular  intervals  as  shown  in  Fig.  4.  This  technique  was  also 
used  for  organic  sorbate  mixtures  because  uv  spectrophotometrlc  detection 
Was  unsuitable  for  multiple  solute  effluents.  After  collection,  a  known 
volume  of  hexane  containing  an  Internal  standard  was  injected  through  the 
bottle  septum  cap  and  a  single  extraction  was  done.  A  portion  of  the 
hydrocarbon  phase  was  removed  by  a  syringe  needle  and  injected  in  a  GC  for 
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Fig.  2.  Schematic  of  chromatography  apparatus 
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Fig.  4.  Manual  collection  of  column  fractions 


Ol 

(j 

<  5 

o  3 

CM 


Hi  (/)  tu 


5  3 
0-1 
CC  O 
u.  u 


analysis.  It  was  later  noted  that  a  similar  method  was  used  by  Patrick  et 
al  (1985)  in  analysis  of  groundwater  samples. 

Selection  of  Flow  Velocity. 

Dispersion  in  LC  columns  causes  peak  broadening  and  is  avoided  in 
analytical  applications  by  using  fast  flow  rates  and  small  particles  with 
thin  layer  coatings  so  local  equilibrium  is  quickly  approached.  This 
involves  Careful  design  of  packing  materials.  It  is  impossible  to  optimize 
flow  conditions  for  aquifer  material  column  packings,  because  many  different 
incompletely  known  particle  types  and  sizes  are  present.  However,  some 
generalizations  remain  applicable  for  aquifer  material  packing.  High  flow 
velocities  give  less  dispersion  and  less  time  for  local  equilibration.  Very 
low  flow  velocities  may  attain  local  equilibrium  but  with  dispersion  greatly 
overriding  sorption  effects.  This  is  apparent  from  the  column  transport 
equation  of  Bastian  and  Lapidus  (1956).  Flow  velocities  chosen  here  (0.001 
to  0.01  cm/sec)  were  in  the  range  for  sedimentary  aquifer  material  given  in 
Todd  (1976).  Practical  column  methods  of  isotherm  measurement  require  that 
dispersion  not  be  too  great  and  that  local  equilibrium  be  attained.  The 
amount  of  dispersion  can  be  estimated  from  the  observed  breakthrough  curves, 
While  proof  of  local  sorption  equilibrium  is  difficult.  If  there  were  no 
other  difficulties  with  column  measurements,  this  proof  could  be  obtained 
from  demonstrated  identity  of  isotherms  measured  by  the  batch  and  column 
methods.  Invariance  of  the  sorption  isotherms  calculated  by  the  column 
methods  with  flow  velocity  is  a  necessary  but  not  sufficient  for  local 
equilibrium.  The  "correctness"  of  column  dynamic  isotherms  can  only  be 
verified  by  comparison  with  batch  isotherm  values,  and  literature  indicates 


column  and  batch  methods  may  not  agree  exactly  (e.g.,  Bllkert  and  Rao, 

1985). 

Flow  velocities  were  too  slow  to  use  the  frontal  analysis  by 
characteristic  point  (FACP)  method  (Jacobson  et  ax,  198^4)  for  isotherm 
measurement.  Breakthrough  curves  from  the  present  study  indicated 
considerable  dispersion,  and  the  FACP  method  does  not  contain  a  means  of 
correcting  calculated  isotherms  for  dispersion.  An  empirical  FACP 
correction  has  been  proposed  by  Zverev  et  al  (197^),  but  has  not  yet  been 
supported  by  mathematical  analysis  of  the  column  transport  process.  There 
may  be  considerable  departure  from  local  equilibrium  at  high  flows,  but  this 
could  not  be  checked  because  calculated  sorption  coefficients  by  the  frontal 
analysis  method  were  found  independent  of  flow  rate  below  that  producing 
column  retainer  blockage.  Blockage  was  indicated  by  a  sudden  rise  of  back 
pressure  from  less  than  100  psi  to  the  column  structural  safety  limit  (2700 
psi)  specified  by  the  manufacturer. 

Measurement  porosity  and  system  dead  volume. 

Frontal  analysis  requires  measurement  of  the  column  and  system  dead 
volume  which  is  the  sum  of  the  column  void  volume  and  the  volume  of  tubing, 
photometric  cells,  column  ends  and  fittings.  The  latter  are  known,  so 
the  task  of  measurement  of  dead  volume  becomes  that  of  determining  the 
column  void  volume,  which  is  also  used  in  calculation  of  porosity.  The 
Obvious  method  is  to  weigh  the  dry  packed  column  empty  and  again  full  of 
solvent.  Weight  difference  and  the  solvent  density  give  the  void  volume 
from  which  one  can  calculate  porosity.  Difficulties  arise  if  the  solvent 
phase  is  compressible  and  applied  pressure  is  above  atmospheric.  A 
Completely  dry  packing  cannot  be  obtained  without  subjecting  clay  minerals 
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and  organic  substances  to  temperatures  where  they  are  irreversibly  altered. 
Packing  density  may  vary  after  solvent  is  passed  and  pressure  is  cycled 
several  times.  Irregular  and  multisized  particles  rearrange  and  sort. 
Aquifer  materials  used  here  were  air  dried  for  24  hours  at  60°C,  and  columns 
were  dry  filled  with  vertically  tapping  as  suggested  by  Snyder  and  Kirkland 
(1979).  Determination  of  the  column  void  volume  is  a  topic  of  active 
research.  Some  of  the  problems  are  reflected  by  the  work  of  McCormick  and 
Karger  (1980)  and  Knox  and  Kalizan  (1985).  The  frontal  analysis  method  adds 
a  constraint  on  dead  volume  determination.  Column  void  volume  by  the 
weighing  method  should  agree  exactly  with  that  calculated  from  the  column 
flow  rate  and  the  observed  breakthrough  curve  of  an  unretained  solute. 

A  mechanical  problem  with  column  packing  becomes  apparent  when  one 
tries  to  assure  that  there  is  no  dead  space  between  the  packing  retainers 
and  the  packing,  or  that  the  packing  bulk  density  is  homogeneously 
distributed  in  the  column  space.  Techniques  in  practical  liquid 
Chromatography  have  been  developed  to  deal  with  this.  Commercial  units  can 
compress  plastic  LC  columns  after  they  are  packed.  Some  columns  have 
movable  end  pieces  that  can  be  adjusted  to  conform  to  and  slightly  compress 
the  packing  bed.  Neither  method  is  applicable  here  because  organic  polymers 
must  be  used  in  their  construction.  A  satisfactory  solution  to  this  problem 
has  not  yet  been  found  using  acceptable  construction  materials,  i.e. 
stainless  steel  and  glass. 

Selection  of  an  unretained  solute  for  frontal  analysis. 

Several  substances  have  been  considered  as  unretained  solutes  for 
groundwater  contaminant  transport  research;  Br  or  Cl  ionic  salts  with 
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alkali  and  alkaline  earth  cations;  benzoic  acid  in  anion  form, 

(tritiated  water);  (deuterated  water).  Freon  12  and  Freon  23. 

Betson  et  al  (1985)  describe  the  use  of  several  of  these  tracers  in  column 
and  field  studies  of  groundwater  transport.  The  halides,  Br  and  Cl  are 
most  often  used  on  the  assumption  that  small  anionic  species  will  not  be 
retained  by  the  negative  charges  on  aquifer  material  surfaces  at  normal  pH. 
Data  given  in  Nkedi-Kizza  et  al  (1983)  for  soil  column  breakthrough  of 
pulses  containing  ^H^O,  '*^Ca  and  **C1  labels  in  CaCl2  solutions  indicated 
that  Ca^^  ions  are  more  retained  than  ’®C1  or  ’H^O,  which  are  themselves 
nearly  similarly  retained.  Thus,  refractive  index  measurement  is  not  a 
perfect  breakthrough  indication  for  halide  salt  solutions  because  the  cation 
is  retained,  probably  in  a  complex  ion  exchange  reaction  with  H  and  other 
cations  on  the  aquifer  material.  Chloride  analysis  of  the  column  effluent 
gives  a  better  approximation  to  an  unretained  solute.  The  ideal  non- 
retained  solute  should  not  interact  with  sorbent  surfaces.  ’HjO  protons 
will  exchange  with  Si-OH  protons  and  organic  functional  group  protons,  and 
may  thus  be  retained.  The  chemistry  of  the  aquifer  material  sorbent  is 
altered  by  the  presence  of  a  dissolved  salt  in  a  complex  and  unknown  way,  so 
halide  salts  or  other  ions  are  not  ideal  unretained  substances. 

Justification  for  the  use  of  Cl  here  is  based  on  the  above-noted  conclusion 
of  similarity  to  transport.  Nkedi-Kizza  et  al  (1982)  have  shown  that 

’®C1  is  very  slightly  retained  relative  to  tritiated  water  on  soils  at  low 
flow  rates,  but  this  difference  was  neglected  here.  Tritiated  water  could 
not  be  used  at  the  time  due  to  license  limitations.  Deuterated  water  is 
probably  suitable,  but  expensive  to  analyze.  The  noble  gases  and  freons 
require  expensive  analysis  and  special  sample  handling  methods  to  avoid 
volatility  losses. 
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Proof  of  p.on-retention  must  await  greater  knowledge  of  the  sorbent 
surface  and  its  chemical  interactions  with  these  solutes.  Nkedi-Klzza  et  al 
(1982)  allude  to  this  problem  in  their  consideration  of  pH  interdependence 
with  salt  ion  sorption  for  CaCl^,  but  of  course  cannot  explain  behavior  in 
this  terribly  complex  situation. 

Imposition  of  a  step  increase  in  unretained  solute  concentration  must 
alter  the  surface  of  the  sorbate  in  ways  which  are  now  undeterminable. 

There  is  no  ideal  unretained  reference  solute.  If  the  solute  chosen  is  not 
perfectly  unretained,  dead  volume  calculated  will  be  too  high  and  the 
sorption  coefficients  measured  by  frontal  analysis  will  be  underestimated. 

It  appears  that  is  the  best  available  "unretained"  substance  for 

aquifer  material,  but  this  is  not  confirmable  at  present. 

Column  design  requirements. 

The  requirements  of  columns  to  be  used  with  aquifer  material  packings 
are  stringent: 

a.  Dead  space  at  column  ends  must  be  negligible  relative  to  packing 
void  volume. 

b.  Column  packing  retainer  must  retain  all  particles  of  packing,  but 
not  block  flow. 

c.  Column  and  system  materials  must  not  sorb  solutes  to  be 
investigated . 

d.  Columns  must  be  free  from  corrosion  by  salt  solution  in  the  presence 
of  some  oxygen,  as  corrosion  products  of  metals  are  good  sorbents. 

e.  Columns  must  withstand  applied  pressure  to  at  least  100-200  psi . 

f.  Wall  and  end  effects  must  be  minimized  by  column  design. 

Heqijirements  c,  d,  and  e  apply  to  chromatographic  system  tubing, 
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connections,  and  measurement  cells.  The  all  stainless  steel  system  used 
here  does  not  meet  requirements  b  and  d.  An  ideal  system  meeting  these 
requirements  for  solute  listed  would  be  all  heavy  borosilicate  glass,  with 
sufficiently  fine  glass  frit  packing  retainers,  glass  connections,  and  glass 
capillary  tubing.  It  should  have  no  plastic  seals  or  organic  lubricants. 

No  such  chromatographic  system  is  commercially  available,  but  designs  are 
being  developed  at  this  laboratory.  It  is,  however,  relatively  simple  to 
design  an  ail-glass  ampoule  system  for  performing  batch  equilibration 
sorption  measurements  on  systems  of  aqueous  solutions  and  aquifer  materials. 

Consideration  of  Relative  Simplicity  for  Column  and  Batch  Isotherms. 

Batch  isotherms  require  constant  temperature  and  measurement  of  liquid 
and  solid  phase  concentrations  and  amounts.  Column  isotherms  require 
Constant  temperature  and  flow  rate,  along  with  measurement  of  packing 
porosity  and  of  two  solute  breakthrough  concentrations  as  functions  of  time. 
Column  isotherm  measurement  is  more  complex  and  this  implies  more  potential 
sources  of  error  in  sorption  coefficients.  It  is  not  satisfying  to 
use  complex  methods  and  costly  equipment  to  obtain  data  that  can  be 
collected  more  simply.  Accordingly,  batch  techniques  appear  preferable  to 
column  techniques  for  sorption  measurement.  A  batch  system  is  better 
definable  than  a  dynamic  flowing  column  system  because  knowledge  of 
hydrodynamics  of  the  former  system  is  not  necessary  for  batch  system 
thermodynamic  description. 

Competitive  sorption. 

Competitive  sorption  was  considered  in  multiple  sorbate  column 
experiments  done  with  1  methylnaphthalene-1 ,2  di chlorobenzene  solutions,  but 
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as  for  single  sorbates,  one  point  K  values  were  obtained  due  to  column 
measurement  time  limitations.  The  time  required  for  one  complete  isotherm 
with  seven  points  by  frontal  analysis  is  about  thirty  days  under  the 
conditions  used.  Complete  batch  isotherms  on  these  materials  can  be  done 
much  more  rapidly  (in  about  two  days)  and  were  done  on  each  of  the  aquifer 
materials  with  1,2  dichlorobenzene  and  with  1  raethylnaphthalene ,  and  on 
Carswell  aquifer  material  with  1  methylnaphthalene-1 ,2  di chlorobenzene 
mixtures.  All  batch  isotherms  were  linear  at  concentrations  used  in  the 
column  experiments,  so  linear  sorption  coefficients  are  presented. 

Competitive  sorption  is  defined  by  Helfferich  and  Klein  (1970)  as  an 
interference  resulting  from  competition  of  the  different  sorbate  molecules 
in  solution  for  the  limited  capacity  (or  number  of  sites)  on  the  sorbent. 
Addition  of  different  molecules  to  the  mobile  phase  of  fixed  composition  in 
another  solute  will  reduce  its  solvent  phase  concentration  due  to 
competition.  Competitive  sorption  is  the  basis  for  displacement 
Chromatography,  and  Jacobson  et  al  (1987)  have  presented  an  excellent  work 
on  measurement  of  competitive  adsorption  isotherms  by  frontal  chromatography 
that  is  relevant  to  groundwater  transport  of  contaminants.  Competitive 
sorption  of  components  of  a  sorbate  mixture  can  be  identified  by  changes  in 
the  isotherm  of  a  given  component  caused  by  the  presence  of  other  sorbates. 
Jacobson  et  al .  (1987)  showed  competitive  sorption  of  phenol  and  p-cresol 
from  water  on  octadecyls^l ica ,  but  this  competitive  sorption  did  not  occur 
below  5  milimole/liter ,  At  the  low  concentrations  presented  in  Table  2, 
isotherms  for  mixtures  did  not  show  competition  and  were  linear  as  noted 
above . 

Competitive  sorption  had  not  been  demonstrated  in  literature  for 
ground  water  contaminant  transport,  and  was  not  detected  by  batch  or  column 
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experiments  done  here.  Absence  of  competition  is  explained  by  the  low 
organic  solute  concentrations  in  contaminated  groundwater  and  by  the 
predominance  of  non-specific  interactions  (o.g.,  van  der  Waals  forces,  etc.) 
in  uptake  of  non-polar  sorbates  by  aquifer  materials.  If  the  column 
experiment  results  of  Jacobson  et  al  (1987)  are  applied  to  groundwater 
aquifer  material  situations  by  assuming  octadecylsil ica  "similar"  to  an 
aquifer  material,  it  appears  that  non-linearity  and  sorbate  competition  do 
not  occur  below  phenol  concentrations  of  about  400  mg/liter.  This  is  higher 
than  the  water  solubility  of  any  sorbate  in  this  program,  and  a  factor  of 
ten  higher  than  the  water  solubility  of  complete  JP-4  fuel  (see  for  example. 
Smith  et  al  1981).  Competitive  sorption  may  occur  in  the  vicinity  of 
environmental  releases  of  mixtures  of  large  amounts  of  highly  water  soluble 
polar  organic  compounds,  but  apparently  does  not  occur  for  liquid 
hydrocarbon  fuels  or  chlorinated  solvents. 
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RESULTS 


Results  of  column  sorption  of  QDCB  and  ICHaN  on  Aqulfci’  Materials. 

Conditions  and  results  of  frontal  chromatography  of  ODCB  and  ICHjN 

solutions  on  aquifer  materials  are  presented  in  Table  2.  Sorption 

Coefficients  for  the  same  systems  from  batch  equilibrations  are  presented 

here  (from  Stauffer,  MacIntyre,  and  Wickman  1989)  for  comparison  with  the 

column  derived  sorption  coefficients.  Trends  observed  in  batch  sorption 

coefficients  were  followed  by  the  column  sorption  coefficients,  but  the 

column  coefficients  were  always  lower  than  batch  coefficients.  Means  and 

standard  deviations  of  the  ratio  (K  from  V.)/K  ^  ^  ^  for  each  sorbate 

1  batch 

averaged  over  all  sorbents  are: 

ODCB,  X  =  0.il9  and  a  =  +0.17 
ICHjN,  X  =  0.45  and  a  =  _+0.12. 

This  implies  an  empirical  relation  such  as: 

column  batch 

This  discrepancy  between  methods  is  probably  not  experimental  error. 
Bilkert  and  Rao  (1985)  observed  a  similar  factor  between  column  and  batch 
sorption  coefficients  for  nematicides  on  soils.  They  noted  that  this  may  be 
a  consequence  of  failure  to  reach  local  equilibrium  in  the  column.  Two  step 
kinetics  have  been  observed  in  batch  approach  to  equilibrium  data,  and  it 
may  be  that  the  second  slower  step  is  not  completed  during  the  contact 
Period  in  columns.  In  this  case  the  observed  K  might  be  a  function  of 
linear  flow  rate,  but  linear  flow  velocities  used  here  were  varied  by  a 
factor  of  ten  from  0.01  cm/sec  to  0.001  cm/sec  without  apparent  effect  on  K. 
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Table  2.  Sorption  Chromatography  Results  on  Aquifer  Materials  (from  Stauffer, 
MacIntyre,  and  Wickman  1989). 

Sorbent  Q  d  Sorbate  C  W  e 


Carswell 

3 

6.2 

ODCB 

6 

6.2 

ICH3N 

tf 

3 

6.2 

ICH3N 

tt 

1  .5 

6.2 

ODCB 

tt 

1  .5 

6.2 

ICH3N 

5 

21  .2 

ICH3N 

!? 

5 

21  .2 

ICH3N 

5 

21  .2 

ICH3N 

11 

5 

21  .2 

ICH3N 

Tinker 

3 

6.2 

ICH3N 

tf 

3 

6.2 

ICH3N 

tf 

3 

6.2 

ODCB 

ff 

3 

6.2 

ODCB 

tf 

3 

6.2 

ICH3N 

tf 

3 

6.2 

ODCB 

If 

3 

6.2 

ODCB 

ff 

3 

6.2 

ICH3N 

11 

3 

6.2 

ICH3N 

tf 

3 

6.2 

ICH3N 

ft 

3 

6.2 

ICH3N 

ff 

3 

6.2 

ODCB 

ff 

3 

6.2 

ODCB 

tf 

3 

6.2 

1C'i3N 

Borden 

1  5 

6.2 

IC.’^N 

ft 

1  .5 

6.2 

ICH3N 

ff 

1  .5 

6.2 

ICH3N 

Lula 

0.75 

4 . 6 

1CH,N 

tf 

0.75 

4 . 6 

ICH3N 

Blytheville 

5 

21  .2 

ICHjN 

ff 

5 

21  .2 

ODCB 

tf 

5 

21  .2 

ICH3N 

If 

5 

21  .2 

ODCB 

Barksdale 

10 

21  .2 

ICHjN 

ff 

10 

21  .2 

ODCB 

10 

21  .2 

ODCB 

1/4 

13.48 

0.361 

1  .78 

2.79 

1/4 

13.48 

0.361 

1 .78 

2.79 

1/4 

13.48 

0.361 

1  .78 

2.79 

1/4 

1  3.48 

0.361 

1  .78 

2.79 

1/4 

13.48 

0.361 

1 .78 

2.79 

1/4 

153.24 

0.368 

1  .74 

2.79 

1/4 

153.24 

0.368 

1  .74 

2.79 

1/4 

153.24 

0.368 

1  .74 

2.79 

1/4 

153.24 

0.368 

1  .74 

2.79 

1/4 

13.72 

0.745 

1  .81 

2.77 

1/8 

13.72 

0.345 

1  .81 

2.77 

1/2 

13.72 

0.345 

1  .81 

2.77 

1/2 

13.72 

0.345 

1  .81 

2.77 

1/4 

13.72 

0.345 

1  .81 

2.77 

1/2 

13.72 

0.345 

1.81 

2.77 

1/8 

13.72 

0.345 

1  .81 

2.77 

1/4 

13.72 

0.345 

1 .81 

2.77 

1/16 

13.72 

0.345 

1 .81 

2.77 

1/8 

13.72 

0.345 

1 .81 

2.77 

1/2 

j  • '  ‘ 

■..345 

1 .81 

2.77 

1/2 

1  p.  f2 

0.345 

1  .81 

2.77 

1/4 

13.72 

0.345 

1 .81 

2.77 

1/2 

13.72 

0.345 

1 .81 

2.77 

1  /? 

13.68 

0.326 

1 .81 

2.67 

1/h 

10.68 

0.326 

1 .81 

2.67 

1/8 

13.68 

0.326 

1 .81 

2.67 

1/2 

7.42 

0.336 

1  .79 

2.69 

1/2 

7.42 

0.336 

1 .79 

2.69 

1/4 

152.50 

0.340 

1  .73 

2.62 

1/2 

152.50 

0.340 

1.73 

2.62 

1/4 

152.50 

0.340 

1 .73 

2.62 

1/8 

152.50 

0.340 

1 .73 

2.62 

1/3 

157.75 

0.351 

1 .79 

2.75 

1/4 

157.75 

0.351 

1.79 

2.75 

1/3 

157.75 

0.351 

1 .79 

2.75 

Q  is  column  flow  rate  (ml/h) 
d  is  column  diameter  (m) 

C  is  sorbate  concentration  as  fraction  of  saturation  concentration 
W  is  packing  weight  (g) 
e  is  packing  porosity 
p  is  bulk  density 
p  is  packing  material  density 
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Table  2  (continued) 


Sorbent 

V, 

Vi 

V/ 

k  from 
.  V,' 

< 

0 

3 

k  Batch 

Detection 
for  V ,  ' 

Carswell 

3.91 

2.81 

3.01 

0.07 

0.08 

0.18 

RI 

T! 

:  .bG 

2.81 

3.28 

0.25 

0.29 

0.45 

RI 

IT 

6 . 06 

2.81 

2.93 

0.23 

0.24 

0.45 

RI 

Tf 

3.99 

2.81 

3.10 

0.07 

0.09 

0.18 

RI 

IT 

7.17 

2 . 8l 

3-05 

0.31 

0.32 

0.45 

RI 

»» 

59.8 

32.5 

42.0 

0.12 

0.18 

0.45 

[Cl'] 

ft 

68.3 

32.5 

37.1 

0.20 

0.23 

0.45 

RI 

tt 

67.1 

32.5 

32.0 

0.23 

0.23 

0.45 

[Cl  ] 

73.5 

32.5 

32.0 

0.27 

0.27 

0.45 

[Cl  ] 

Tinker 

5.18 

3.00 

3.02 

0.16 

0.16 

0.44 

RI 

Tt 

5.18 

3.00 

3.02 

0.16 

0.16 

0.44 

RI 

T! 

3.81 

3.00 

3.04 

0.06 

0.06 

0.10 

RI 

IT 

3.88 

3.00 

3.31 

0.04 

0.04 

0,10 

RI 

»1 

4.48 

3.00 

3.09 

0.10 

0.11 

0.44 

RI 

IT 

3.79 

3.00 

3  07 

0.05 

0.06 

0.10 

RI 

11 

3.75 

3.00 

3.03 

0.05 

0.05 

0.10 

RI 

tl 

5.19 

3.00 

3.06 

0.16 

0.16 

0.44 

RI 

TT 

5.42 

3.00 

2.99 

0.18 

0.18 

0.44 

RI 

TT 

5.20 

3.00 

^.93 

0.17 

0.16 

0.44 

RI 

Tl 

4.57 

3.00 

2.93 

0.12 

0.11 

0.44 

RI 

T! 

3.31 

3.00 

2.93 

0.03 

0.02 

0.10 

RI 

TT 

3.75 

3.00 

3.07 

0.05 

0.05 

0.10 

RI 

TT 

5.92 

3.00 

3.^7 

0.18 

0.21 

0.44 

RI 

Borden 

5.75 

2.46 

2.90 

0.21 

0.24 

0.65 

RI 

TT 

5.95 

2,46 

2.91 

0.22 

0.26 

0.65 

RI 

TT 

6.05 

2.46 

2.79 

0.24 

0.26 

0.65 

RI 

Lula 

3.12 

1  .39 

1  .73 

0.19 

0.23 

0.51 

RI 

TT 

2.80 

1.39 

1 .46 

0.18 

0.19 

0.51 

RI 

Blytheville 

182.3 

29.0 

28.75 

1  .00 

1  ...1 

3 . 06 

[Cl'] 

TT 

82.8 

29.0 

30.00 

0.35 

0.35 

1  .28 

[Cl  ] 

TT 

295.0 

30.0 

— 

-- 

1.74 

3.06 

— 

T! 

110.5 

30.0 

— 

-- 

0.53 

1  .28 

— 

Barksdale 

128.46 

31  .b 

36.32 

0.58 

0.62 

0.96 

[Cl'] 

IT 

66 .  bb 

31  .0 

39.26 

0.18 

0.23 

0.30 

[C1_] 

69.57 

31  .0 

40.25 

0.18 

r\i 

0 

0.30 

[Cl  ] 

Vj  is  retention  volume  of  sorbate  (mi) 

'.'i  is  column  dead  volume  (ml)  from  weight  data 

Vj'  is  column  dead  volume  (ml)  from  unretained  substance  breakthrough  curve 
k  is  calculated  sorption  coefficient  (liters/kg)  assuming  linear  isotherm 
RI  indicates  refractive  index  measurement 
[Cl  ]  indicates  chloride  titration 
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Rao  and  Jessup  (1983)  speculated  upon  possible  reasons  that  sorption 
equilibrium  is  not  reached  in  columns  at  low  flow  rates  that  should  yield 
sufficient  contact  time.  They  considered  several  models  involving  liquid 
filled  pore  dead  space,  diffusion  through  organic  coatings  on  sorbent,  and 
multiple  kinds  of  sorption  sites.  No  conclusions  were  reached,  and  the 
present  work  does  not  resolve  the  matter.  It  is  also  possible  that  the 
severe  turbulence  and  energetic  collisions  between  sorbent  particles  in  a 
batch  isotherm  measurement  may  abrade  particles  and  create  new  surface  area 
and  sites  for  rapid  equilibration,  perhaps  to  a  different  state  than  exists 
in  columns.  The  author  has  earlier  noted  several  difficulties  with  liquid 
column  chromatography  measurement  of  sorption  coefficients;  some  of  these 
can  lead  to  apparently  low  K  with  respect  to  batch  K  (e.g.,  particle  loss, 
packing  channeling).  Rao  and  Jessup  (1983)  state  that  further  research  on 
the  disagreement  between  batch  and  column  sorption  is  needed,  and  the 
present  author  concurs. 

The  observation  of  a  general  factor  of  about  0.5  almost  independent  of 
particular  soil,  aquifer  natural  or  sorbate  is  tantalizing.  An  explanation 
might  be  obtained  from  a  set  of  concurrent  batch  and  column  sorption 
measurements  of  selected  slightly  polar  organic  compounds  on  uniform 
particle  size  synthetic  sorbents  of  somewhat  known  composition.  Questions 
to  be  answered  include: 

1)  Does  the  factor  hold  for  an  inorganic  sorbate,  such  as  geothite  or 
kaolinite? 

2)  Does  the  factor  hold  for  a  C^g  bonded  phase  spherical  pellicular 
packing  sorbate? 

3)  Do  functional  group  substitutions  on  0,8  by  -CN,  -NH^  etc  in  the 
bonded  phase  affect  the  factor? 
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4)  Does  the  factor  hold  both  for  pellicular  and  microporous  bonded 
phase  particles? 

5)  Is  the  factor  dependent  on  the  sorbent  particle  size  in  any  of  the 
above  experiments? 

The  failure  of  column  K  data  to  agree  with  batch  K  data  precludes  the 
initially  proposed  thermodynamic  analysis  of  column  data  using  aquifer 
materials.  There  are  no  apparent  advantages  to  column  measurement  of  soil 
and  aquifer  material  sorption,  and  it  is  recommended  that  batch  isotherm 
data  be  used  for  comparison  of  aquifer  sorption  purposes,  at  least  until 
differences  between  batch  and  column  methods  are  understood.  Kinetic 
investigations  of  sorption  by  column  methods  have  been  proposed  by  Rao  and 
Jessup  (1983).  but  seem  fraught  with  complexity  and  do  not  permit  unique 
interpretations  of  results.  Kinetic  measurements  yield  unique  sorption 
mechanisms  only  when  intermediates  are  observed  microscopically  or 
spectroscopically.  The  author  believes  kinetics  of  sorption  in  systems 
studied  must  be  inerpreted  empirically,  and  that  there  is  now  no  best  model 
for  sorption  kinetics  measurement.  Instead  there  are  different,  but 
mathematically  similar  methods  using  approximate  and  incomplete  models  of 
the  sorption  process. 

Table  3  presents  arithmetic  average  sorption  coefficients  for  each 
aquifer  material  by  both  batch  and  column  methods.  ODCB  was  not  run  on 
Borden  and  Lula  materials.  The  order  of  uptake  by  sorbents  is  similar 
within  experimental  error  for  the  two  methods. 
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Table  3.  Sorption  coefficients  for  1CH  N  and  ODCB  on  aquifer  materials  by 
column  and  batch  techniques.  ^ 

COLUMN  BATCH 


ICH^N 

ODCB 

ICH^N 

ODCB 

Tinker 

0.18 

0.05 

0.44 

0.10 

Carswell 

0.25 

0.08 

0.45 

0.18 

Blytheville 

1.38 

0.44 

3.06 

1  .28 

Borden 

0.26 

— 

0.65 

0.30 

Lula 

0.21 

— 

0.51 

0.14 

Barksdale 

0.62 

0.23 

0.96 

0.30 

It  can  be 

argued  that,  if 

reality  and  correctness  of  the  factor  near 

0.5  is  accepted 

,  there  is  good 

agreement 

between  the 

sorption  coefficients 

for  the  aquifer  material  which  were  determined  by  independent  analytical 
methods.  However,  application  of  this  factor  without  understanding  is  not 
satisfying. 

Work  with  sorbate  binary  mixtures  was  limited  due  to  the  need  for  GC 
analysis  of  many  fractions  from  each  column  run.  Mixture  studies  were 
further  restricted  when  discrepancies  were  noted  between  batch  and  column 
methods,  and  time  was  allocated  to  establishing  the  magnitude  of  this 
difference.  The  last  ODCB  and  ICH^N  sorption  coefficients  on  Blytheville 
material  given  in  Table  2  were  measured  on  an  ODCB  -  ICH3N  mixture  such  that 
column  influent  water  was  1/4  saturated  with  ICH3N  and  1/8  saturated  with 
ODCB.  The  sorption  coefficients  obtained  were  slightly  higher  than  those 
Obtained  in  single  solute  uptake  measurements,  but  there  was  insufficient 
data  on  Blytheville  material  to  establish  statistically  significant 
differences . 
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Breakthrough  curves  for  JP-4  solution  on  aquifer  material. 


The  purpose  of  these  experiments  was  to  demonstrate  the  chromatography 
of  JP-4  components  under  emulated  aquifer  conditions.  JP-4  was  equilibrated 
With  water  in  the  volume  ratio  of  1:100  to  produce  a  saturated  solution. 

This  aqueous  solution  was  diluted  by  an  equal  volume  of  water  and  loaded 
into  a  syringe  pump.  Chromatography  was  done  under  conditions  similar  to 
those  used  for  simple  hydrocarbon  solutions;  flow  rates  5  ml/h,  column  25  cm 
X  21.2  mm  ID,  packing  with  Blytheville  aquifer  material,  dead  volume 
calculation  from  Cl  breakthrough  curve  when  CaClj  concentration  was  stepped 
from  0.01  to  0.015M.  Breakthrough  of  components  of  the  solution  was 
measured  by  collecting  fractions  over  a  2-hour  period  in  septum  capped 
bottles  kept  at  0°C  in  an  ice  bath,  extracting  the  content  with  1  ml  hexane 
containing  tri- isopropyl  benzene  internal  standard,  and  injecting  the 
extract  into  a  gas  chromatograph  with  capillary  column  to  analyze  each 
component . 

The  breakthrough  curves  were  generated  from  discrete  points.  Graphical 
integration  was  done  on  manually  smoothed  data,  and  the  results  checked  by 
numerical  integration  in  wnich  areas  were  calculated  for  trapezoids 
corresponding  to  each  data  point. 

Establishment  of  breakthrough  curves  for  JP-4  components  by  this  method 
was  only  possible  for  those  components  that  are  quite  soluble  in  water  and 
abundant  in  the  JP-4.  High  molecular  weight  compounds  (i.e.,  highly 
substituted  benzenes  and  naphthalenes)  and  aliphatic  hydrocarbons  did  not 
meet  these  conditions.  Good  breakthrough  data  was  obtained  for  toluene, 
ethylbenzene,  and  the  xylenes.  The  m  and  p  xylenes  were  not  separated 
by  the  gas  chromatographic  method  used,  so  their  sum  was  recorded,  and  used 
with  the  assumption  that  the  two  isomers  also  behaved  similarly  in  liquid 
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chromatography  on  the  aquifer  material.  This  is  supported  by  the 


breakthrough  curves  are  shown  in  Fig.  5.  All  concentrations  are  normalized 
to  those  in  the  feed  solution.  Calculated  K  values  were: 


K 

Boiling  Point  (°C) 

tol uene 

o 

o 

110.6 

ethylbenzene 

0.39 

136.2 

m+p  xylene 

O.iiS 

m  1 39.1  ,  p  138.1 

o-xylene 

0.46 

144.4 

The  K  values  increase  with  boiling  point  and  with  frontal  elution  time  in 
the  same  order  as  peaks  for  the  compounds  appeared  in  GC  analysis  of 
extracts  of  fractions  containing  all  five  solutes.  This  trend  is  expected 
but  cannot  be  explained  in  terms  of  a  sorption  mechanism  on  the  aquifer 
material.  Identities  of  the  compounds  were  confirmed  by  GC-MS  of  the 
extracts.  Sorption  coefficients  observed  here  are  in  reasonable  accord  with 
those  measured  for  ODCB  and  iCHjN,  0.53  and  1 .74  respectively,  on 
Blytheville  material  under  the  same  conditions.  Results  in  Fig.  6  show 
chromatography  which  might  be  anticipated  for  a  JP-4  release  to  groundwater. 
This  situation  has  not  yet  been  demonstrated  here  in  field  studies.  Patrick 
et  al  (1986)  have  injected  water  solutions  of  simple  aromatic  hydrocarbon 
mixtures  In  tests  at  the  Borden  site,  but  have  not  yet  injected  liquid 
hydrocarbon  mixtures  or  petroleum  products.  Chromatographic  separation  must 
occur  for  complex  mixtures,  as  demonstrated  here  in  laboratory  emulation, 
but  sampling  and  analyses  are  much  more  difficult  than  with  simple  mixtures. 
Burris  and  MacIntyre  (1986)  investigated  solution  of  complex  hydrocarbon 
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mixtures  in  water  and  noted  that  the  approximate  concentration  of  component 
compound  in  water  is  given  by  where  is  the  mole  fraction  of 

i  in  the  hydrocarbon  mixture  and  C°  iS  the  solubility  of  the  pure 
hydrocarbon.  Each  component  i  of  a  fuel  mixture  is  present  at  low  X^,  and 
will  be  small  for  aliphatic  or  high  molecular  weight  aromatic  compounds, 
so  to  be  measured  is  small.  This  means  that  large  ground  water  sample 
volumes  must  be  taken  in  the  field  and  that  laboratory  column 
chromatographic  emulations  of  field  situations  require  large  sample  volumes 
which  are  limited  by  the  practical  size  of  chromatographic  columns.  The 
author  estimates  that  one  liter  samples  should  be  taken  for  field  release 
studies  similar  to  the  Borden  work  using  jet  fuels,  and  that  these  samples 
will  require  analysis  using  both  solvent  extractions  and  purge  and  trap 
isolation,  followed  by  GCMS  and  GC-FID.  Analysis  of  high  molecular  weight 
fuel  components  in  groundwater  has  been  ignored  in  most  reported  field 
studies  and  site  investigations. 

Emulation  and  Field  Study  of  Nonlinear  Sorption  in  Contaminant  Transport  by 
Groundwater . 

Groundwater  contaminant  transport  models  that  are  applied  include 
sorption  by  material  in  the  saturated  zone,  and  are  currently  based  on  the 
assumption  of  a  linear  adsorption  isotherm.  Sorption  is  then  expressed  in 
transport  equations  by  a  retardation  factor,  for  example  as  presented  in  the 
compilation  by  Van  Genuchten  (1982).  Contaminant  distributions  predicted 
from  these  equations  will  be  in  error  when  sorption  isotherms  are  nonlinear, 
and  other  factors  such  as  biodegradation  and  vaporization  may  be  Invoked  to 
explain  the  apparent  differences  between  model  prediction  and  field 
observation.  Nonionic  organic  compounds  may  have  nonlinear  isotherms  on 
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soils  and  aquifer  materials,  and  some  examples  are  shown  in  Mingelgrin  and 
Gerstl  (1983).  Transport  equations  with  linear  retardation  factors  have 
been  used  as  approximations  because  inclusion  of  nonlinear  isotherms 
produces  equations  more  difficult  to  solve  analytically.  A  general  analytic 
solution  of  the  transport  equation  for  any  nonlinear  isotherm  has  not  yet 
been  obtained,  but  considerable  progress  has  recently  been  made  with 
analytic  solutions  in  one  dimension  for  parabolic  and  Langmuir  isotherms. 

This  discussion  applies  some  results  from  nonlinear  chromatography  to 
groundwater  contaminant  transport  and  proposes  field  tests  to  investigate 
the  importance  of  isotherm  nonlinearity  in  determining  contaminant 
distributions  in  aquifers.  Field  studies  of  nonlinear  sorption  process?'’  at 
sites  similar  to  the  current  Borden  site  are  suggested,  because  aquifer 
material  in  situ  chemical  environment,  fluid  dynamics  and  sedimentology 
cannot  be  duplicated  in  the  laboratory.  Sposito  et  al  (1986)  noted  the  need 
for  field  experiment  tests  of  contaminant  transport  models  in  their 
discussion  of  stochastic  transport  models.  In  the  following,  some  recent 
results  from  nonlinear  adsorption  chromatography  literature  are  presented 
and  their  relevancy  to  groundwater  contaminant  transport  calculations  is 
considered.  A  field  site  research  program  that  involves  introduction  of 
organic  compounds  having  linear  and  nonlinear  isotherms  on  the  aquifer 
material  is  proposed.  This  will  demonstrate  the  extent  to  which  contaminant 
sorption  isotherm  nonlinearity  produces  contaminant  distributions  and 
transport  velocities  of  transport  that  differ  from  predictions  assuming 
linear  sorption. 

Recent  chromatographic  research  on  nonlinear  sorption  liquid 
chromatography  has  treated  advective-dispersi ve  transport  with  sorption  for 
different  classes  of  isotherms.  Jaulmes  et  al  (1986)  treated  the  case  of 
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second-order  polynomial  isotherms  that  can  closely  approximate  Langmuir  and 
some  other  observed  isotherms.  Arinshtein  (1983-  1986)  has  given  analytic 
solutions  for  advecti ve-dispersive  transport  with  a  Langmuir  isotherm, 
generalizing  initial  and  boundary  conditions  in  the  latter  paper.  Zolotarev 
'1985)  has  given  solutions  for  advecti ve-dispersive  transport  with  S-shaped 
isotherms  for  step  function  inputs  (frontal  chromatography).  This  area  of 
research  is  active  and  more  solutions  for  other  isotherms  and  initial  and 
boundary  condition  are  anticipated  (C.  Horvath,  personal  communication). 

The  present  effort  to  relate  this  fundamental  work  on  sorption 
chromatography  to  groundwater  contaminant  transport  selects  the  formulations 
of  Jaulmes  et  al  (1986)  for  further  consideration,  because  of  the  general 
utility  of  isotherm  fitting  by  polynomials,  and  because  this  work  is 
formulated  for  liquid  chromatography.  Liquid  chromatography  terms  are  the 
same  as,  or  directly  translate  to  those  used  in  the  groundwater  literature. 
When  a  pulse  (square  wave)  input  is  applied  to  a  finite  length  column,  the 
expression : 
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gives  the  concentration  of  the  sorbate  at  the  outlet  as  a  function  of  time. 


Where  U  -  limit  band  progation  velocity  for  zero  sample  size 
t  =  time 


;  =  Z-'Jt,  Z  --  length  along  column 

=  concentration  of  sorbate  in  liquid  phase 
aU"X 

u  =  _ 

2D' 
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a  =  peak  area  parameter 


A 


I  ? 


t 


2k, 


t 


t  T  I 

and  k^  and  k^  are  constants  of  the  parabolic  isotherm 

=  partial  molar  volume  of  sorbate  at  zero  concentration 
D  dispersion  coefficient 


D’  =  D/(l  +  k^) 


All  these  parameters  can,  in  principle,  be  measured  independently  and 
used  to  define  the  output  peak.  Jaulmes  et  al  did  not  do  this,  but  used 
experimentally  determined  peak  shape  and  retention  times  to  determine  a,  p, 
A,  and  D.  This  fitting  procedure  is  conceptually  similar  to  that  used  in 
calibratives  groundwater  transport  models. 

Some  features  -  ohis  solution  are  relevant  to  groundwater  situations. 
Peaks  are  3ymm»  .^rical  for  linear  sorption  isotherms.  For  nonlinear  convex 
(Langmuir-Type)  isotherms,  peak  fronts  are  sharp  and  tails  diffuse.  The 
converse  holds  for  concave  isotherms.  The  center  of  sorbate  mass  of  a  peak 
with  linear  sorption  isotherm  is  independent  of  sorbate  concentration,  but 
is  retarded  for  a  convex  isotherm  and  advanced  for  a  concave  isotherm. 
Groundwater  contaminant  plume  shape  and  the  velocity  of  center  of 
contaminant  mass  thus  depends  upon  the  sorption  isotherm  for  the  contaminant 
on  the  aquifer  material. 

If  we  consider  a  pulse  injection  as  has  been  described  at  the  Borden 
site  by  Mackay  et  al  (1986),  plan  views  of  instantaneous  sorbate 
concentration  like  their  Figure  8  will  be  nearly  elliptical  for  sorbates 
with  linear  isotherms.  If  isotherms  are  nonlinear  and  convex,  a  common 
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situation,  the  sorbate  distributions  for  similar  groundwater  conditions  will 
become  tear-drop  shaped  with  the  elongated  portion  of  the  distribution 
directed  back  up  the  hydraulic  gradient.  The  frontal  boundary  of  the  tear 
drop  will  be  very  sharp  with  the  sharpness  limited  by  dispersion.  The 
Center  of  mass  will  arrive  at  a  given  distance  from  the  source  more  quickly 
than  would  have  been  the  case  if  the  isotherm  were  linear. 

Jaulmes  et  al  (1986)  state  this  as: 


D't. 


t  tj-/  1 

m  —  R 


XC 


m 


where  C  =  maximum  peak  concentration 
m 

t  =  limit  of  retention  time  for  zero  sample  size 
R 

t^  =  retention  time  corresponding  to  maximum  peak  concentration 
L  =  column  length 
X  and  D’  have  been  given  above 


or  as  a  linear  approximation  when  XCm  and  are  small  with  respect  to 

unity: 
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If  the  isotherm  is  convex,  k^,  the  first  derivative  of  the  isotherm 

t 

slope  must  be  less  than  one  so  (1  +  k^)  must  be  positive.  UlCjjj  <<  1  for 
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not  too  high  concentrations,  so  t  is  less  than  one.  That  is,  t  <  t„. 

ra  m  R 


These  assumptions  can  be  confirmed  from  the  data  tabulated  by  the  authors. 

Tne  approximate  solution  for  the  convex  polynomial  isotherm  case 
demonstrated  by  Jaulmes  et  al  is  similar  to  the  exact  solution  for  a 
Langmuir  isotherm  by  Arinshtein  (1983). 

The  theory  of  nonlinear  liquid  chromatography  is  not  completely  worked 
out,  but  is  sufficiently  developed  for  use  in  groundwater  contaminant 
transport  calculations.  Its  use  in  transport  models  requires  determination 
of  isotherms  over  the  concentration  range  occurring  in  an  aquifer.  If  the 
isotherms  are  not  exactly  linear,  the  transport  equation  with  a  constant 
retardation  factor  cannot  validly  be  used. 

A  site  at  Columbus,  Mississippi,  is  presently  under  consideration.  The 
research  proposed  here  is  similar  in  concept  to  the  recent  program  of  Mackay 
et  al  (1986)  on  light  chlorinated  hydrocarbons  and  their  more  recent  work  on 
benzene-toluene-xylene  mixtures  (Patrick  et  al  1986).  These  organic  solutes 
all  had  linear  sorption  isotherms  on  the  aquifer  material.  It  is  proposed 
that  similar  releases  of  two  component  mixtures  of  organic  compounds  be 
made,  with  one  of  the  components  linearly  sorbed  and  the  other  nonlinearly 
sorbed.  The  components  must  be  selected  such  that  they  do  not  degrade 
appreciably  during  the  experiment.  Sampling  wells  and  procedures  should  be 
similar  to  those  now  in  use  with  the  Borden  mixture  studies. 

The  linearly  sorbed  compound  should  be  an  aromatic  hydrocarbon  such  as 
ethylbenzene,  i-propylbenzene  or  1-methylnaphthalene.  Several  candidate 
compounds  have  been  found  to  give  nonlinear  isotherms  on  soils.  Selection 
among  those  that  may  give  Freundlich  or  Langmuir  type  isotherms  on  aquifer 
materials  will  involve  consideration  of  a  number  of  polar  nonionic 
pesticides,  herbicides  and  industrial  chemicals.  The  choice  is  based  on 
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extent  of  departure  from  linearity  of  the  sorption  isotherm  and  on 
resistance  to  biodegradation.  Both  requirements  must  be  met  on  the  site 
aquifer  material.  A  high  concentration  of  the  sorbate  will  be  used  to 
assure  that  effects  of  nonlinearity  of  sorption  on  contaminant  transport  are 
observable,  ( hexachlorocyclohexane )  is  a  candidate  compound  with 

isomers  that  Mingelgrin  and  Gerstel  (1983)  have  reported  to  sorb 
nonlinearly,  and  that  is  somewhat  resistant  to  biodegradation.  The  optimum 
nonlinearly  sorbed  compound  will  be  selected  from  results  of  a  laboratory 
program  in  which  batch  sorption  isotherms  and  biodegradability  of  a  number 
of  slightly  soluble  organic  compounds  will  be  determined  using  the  site 
aquifer  material.  Experiments  will  be  conducted  on  columns  packed  with  site 
aquifer  material,  using  the  sorbates  selected  above.  The  results  will 
indicate  chromatographic  processes  that  should  occur  in  field  tests. 

The  objective  is  to  examine  the  proposition  that  nonlinear  sorption 
behavior  of  a  contaminant  solute  in  groundwater  causes  its  transport  rate 
and  concentration  distribution  to  differ  appreciably  from  that  of  a  similar 
solute  exhibiting  a  linear  isothe. m.  If  the  proposition  is  supported  by  the 
proposed  field  work,  contaminant  transport  model  calculation  should  be  based 
on  a  mass  balance  equation  similar  to  that  given  by  oaulmes  et  al  (1986) 
namely , 
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where  Ug  is  the  velocity  of  the  aquifer  flow.  All  other  terms  in  this 
equation  have  been  defined  above.  One  notes  the  differences  between  this 
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equation  and  a  similar  expression  for  a  linear  isotherm  given  for  example  in 
Bear  1979  (using  the  term  definitions  above)  as 


where  k  is  the  linear  sorption  coefficient  times  c/pb.  The  differences  are 
an  expression  of  isotherm  nonlinearity.  If  isotherms  are  measured  carefully 

t  T  t 

to  obtain  good  ana  values,  the  Jaulmes  et  al  (1986)  equation  can  be 
solved  numerically  for  groundwater  transport  situations  for  given  initial 
and  boundary  conditions.  The  field  work  proposed  will  provide  data  needed 
to  demonstrate  solutions  in  an  unconfined  aquifer. 

The  site  at  Columbus,  Mississippi,  is  presently  instrumented  and 
operated  by  The  Tennessee  Valley  Authority.  Initial  visits  indicate  its 
suitability  of  the  site  and  the  potential  to  begin  field  work  in  Ff  1989. 

If  administrative  actions  are  successful,  collection  of  aquifer  material  for 
laboratory  work  will  be  done  in  FY  89.  A  central  theme  of  this  work  is  to 
Consider  transport  of  organics  in  solution,  alone  and  in  mixtures,  as  a  two 
or  three  dimensional  chromatography  problem.  Since  there  is  relatively 
little  vertical  transport  observed  in  most  groundwater  test  site  studies, 
there  is  a  strjng  analogy  between  transport  from  a  ruptured  fuel  tank  or  a 
dump  site  and  the  transport  and  separation  of  chemicals  on  a  thin  layer 
chromatography  plate  (Helling  1971).  The  chromatographic  literature  may  be 
of  considerable  use  in  improving  understanding  of  groundwater  contaminant 
transport. 

Some  recent  numerical  models  have  been  presented  for  groundwater 
contaminant  transport  with  non-linear  sorption,  and  the  SUTRA  model  of  Voss 
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(1984)  has  been  selected  for  initial  use  in  this  research.  SUTRA  is  a  two 
dimensional  finite  element  model  with  vertical  averaging  and  is  suited  for 
contaminant  transport  evaluation  in  the  saturated  zone  using  data  from  the 
proposed  field  test  program. 
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